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A new model is proposed, in the context of Axelrod's model for the study of cultural dissemination, 
to include and external vector field (VF) which describes the effects of mass media on social systems. 
The VF acts over the whole system and it is characterized by two parameters: a non-null overlap 
with each agent in the society and a confidence value of its information. Beyond a threshold value 
of the confidence there is induced monocultural globalization of the system lined up with the VF. 
Below this value, the multicultural states are unstable and certain homogenization of the system is 
obtained in opposite line up according to that we have called negative publicity effect. Three regimes 
of behavior for the spread process of the VF information as a function of time are reported. 

PACS numbers: 89.75.Fb,87.23.Ge,05.50.+q 



I. INTRODUCTION. 

Agent-Based Models (ABMs) 0, [|| are computer simulations of the local interactions of the members of a population 
which could be plants and animals in ecosystems vehicles in traffic, people in society etc. Locally interactions at 
lower-level give rise to the spontaneously emergence of higher-level organizations whose properties are not possessed 
by the individuals neither directly determined by them. Complex and non-linear phenomena have attracted the 
attention of the scientific community to study the interplay between the lower and higher levels of organizations [f| . 
These models typically consist of an environment or framework in which the interactions occur among some number 
of individuals defined in terms of their behaviors (procedural rules) allowing the tracking of the characteristics of each 
individual through time. 

There are lots of applications to model different aspects of dynamics in society @, 0, HI ■ Specifically, the Axelrod 
model 0, [T3| is an ABM designed to investigate the dissemination of culture among interacting agents in a society. 
In this model, society is represented by a lattice composed by a 2-dimensional array of vectors (agents) with a number 
of entries called "features" . The definition of its culture is given by the set of traits an agent has in its features. The 
Axelrod's model has been exhaustively implemented to study a great variety of problems: the nonequilibrium phase 
transition between monocultural and multicultural states [ll| , the cultural drift driven by noise [l2l Il3j ] , nominal and 
metric features [13, Il5|. p ropaganda flU ] , time evolution dynamics [I?} , the resistance of a society to the spread of a 
foreign cultural traits [18| |. finite size effects the impact of the evolution of the network structure with cultural 
interaction [2(| among others. 

Some works have been done including in the Axelrod model an extra agent acting as a vector field (VF) over the 
whole society with the purpose of simulating a mass media effects [2l|, [52, HH, HiJ ■ In all of them, the interaction 
between the external VF and agents is similar to that between an agent and its neighbors: they interact only if they 
have at least one common trait in their corresponding features. In this formalism, the inclusion of an external field, 
which does not change its values on time, introduces an asymmetry on the lattice, which can now be described as 
composed by two groups of agents: group A where agents have trait (s) in common with the VF and group B whose 
do not. All the interactions can now be classified as follows: agents from group A with VF (A-VF), agents from group 
B with VF (B-VF), between agents from group A (A- A), between agents from group B (B-B) and finally between 
agents from group A and B (A-B). See Fig. [T] and lines therein. The (B-VF) interaction, showed as dash-dot-dot 
line in the figure, is a null interaction because agents from group B do not share traits with the VF. In this way, the 
only opportunity of agent B to acquire one VF trait is through a diffusion mechanism with the combined interactions 
(A-VF) plus (A-B) as pointed out in Fig. [T]with a thick line. In those models it is also included the strength of the 
VF as a probability P of interactions A-VF and B-VF while the probability of interactions A-A, B-B and A-B are 
given then by 1 — P. Therefore the diffusion mechanism, showed with a thick lines in Fig. [TJ has very low probability 
as P increases and agents from group B are set apart from the VF information. Furthermore, the mechanism B-VF 
can be very active (time consuming) but with null effects, and the internal relaxing mechanism B-B is not able to 
drive agents to the final state in an efficient way. Thus, the final absorbing states obtained in those previous models 
consisting in multicultural states when the VF strength is increased is then not surprising [U [2^, [H, [24j • 
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FIG. 1: Former models representation of agents and dynamics interactions in the lattice when the VF is included (represented 
by a big square). Agents from group A are indicated with dots while agents in group B are indicated with crosses. Different 
interactions are represented with different lines. The diffusion mechanism from the VF to agents in group B is represented by 
thick lines. Null direct interaction between VF and agents from group B is represented by dash-dot-dot line. 

Our current interest is to develop a new model for the inclusion of an extra agent acting as a vector field (VF) 
over the whole society to overcome the difficulties achieved by the previous models described above. As mentioned in 
Ref. [2l|, the media information is socially processed through personal networks. Then, models will be more realistic 
is they allow a strong interaction with the VF without loss of interchanges between agents on the lattice. It is also 
important to say that mass media designs its publicity in a clever way. As mentioned by the anthropologist Gregory 
Bateson: to produce a change it is necessary to be different but, at the same time, it is necessary to be "close enough" 
to be taken into account [25j |. When acting over the society, mass media always try to have something in common 
with the people chosen as target of publicity or propaganda. It is designed to offers attractive materials for the whole 
society: news, sports, soap operas, movies, cartoons, music, arts, etc. 

Our goal here is to develop a model to include this effect considering an additional non-zero probability of all agents 
to copy a trait from the VF, even if they do not share any trait of their features. Section[TT]is devoted to that purpose. 
In section [TTTI it is exposed some numerical calculations and finally some conclusions arc outlined in section ITVl 

II. THE MODEL. 

The system consists of L 2 agents as the sites of a square lattice. The state of an agent i is defined as a vector of F 
nominal components called features given by <7; = (u^ , 0£, , o~i F ) which characterize the nominal _F-dimensional 
culture of the corresponding agent. This way, each agent has four nearest neighbors but as the fifth it is introduced 
the VF A4 with nominal features o~m = ...,o~Mt > ■■■i (j Mp)- The VF intents to simulate and external mass 

media or publicity which acts over the whole society. Then, each agent can interact with five agents: its four nearest 
neighbors and the VF M, all with equal probability 1/5. Additionally, each feature <Ji f and (JMf can take any of the 
values in the set {0, l,...,q— 1} which are the corresponding cultural traits of an agent i or the VF M.. Initially, the 
values of the vectors o~i and o~m are randomly and independently set with one of the q F state vectors with uniform 
probability. 

The interaction between different agents is possible only when the two vector have an overlap < / < 1 where the 
overlap between two agent i and j is the number of shared traits and it is given by — S/=i ^ a i s - a i s ■ Here <5 is 

the Kronecker symbol. The probability, which we call here nominal probability, of the interaction between two agents 
is given by p(i,j) — l(i,j)/F. In general, the situation p(i,j) = is possible when the overlap between two agents 
is zero, but the case where the probability between an agent and the VF is zero is not an acceptable situation for 
a publicity (or mass media) which intentionally designs its interaction in such a way that always there are features 
which have traits in common with the agent subject of the influence to guarantee that the connection is active. 

In order to include this important effect in our model, we included some effective features that the VF always 
shares with each agent when interacting besides nominal features with the purpose to simulate phenomenologically in 
a simple way the almost omniscient force of today's publicity or propaganda in mass media that offers something for 
all tastes and ages (magazines, radio programs, TV series, etc). Note that the effective features are related with the 
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dynamics between the agents and the VF while the nominal features are related with the dynamics between agents 
inside the lattice. The specific nature in real society of the effective features is not of importance here. It will be 
different for different agents, but the intention is to take into account the specific design of the publicity that mass 
media does to attract everyone. For e/F < 1 where e is defined as the "effective feature", the VF and the agent share 
more nominal than effective features and the VF constitutes a "perturbation" to the internal interaction between 
different agents in the society. The case e/F > 1 means that there are more effective features that certainly share 
each agent with the VF than the number of nominal features each agent has. Then it is now the society which can 
be considered as a "perturbation" with respect to the more robust dynamics between the VF and the agents. In our 
model, the parameter e not only takes natural values, but also fractional values, as it will be seen later. 

Therefore, the probability of interaction between the external vector and an agent, which we call here extended 
probability, is written as 



or, in dimensionless parameters, 



P(hM) = ±-= r + (I) 



where Ai) is the overlap of the nominal features between agent i and the VF. 

In Fig. [2] are shown the values of the probability p(i, M) as a function of e/F for fourth nominal features (F = 4) 
and different values of l(i,M). The values of l/F = 0.00, 0.25, 0.50, 0.75 and 1. 00 are obtained when the agent 
shares with the vector field 0, 1, 2, 3 and 4 nominal features. As seen, the probability is zero only when there are not 
effective features (e = 0) and the overlap between the agent i and the VF is zero. In contrast, in all the other cases 
the probability p(i, Ai) is always different from zero. For e = the values for the case with no effective features are 
recovered. As expected, the probability is lager for larger values of the effective features e for a given value of Ai) 
and also increases for larger values of the overlap Ai) at a given number of e. Finally, when the agent i and the 
VF share all the nominal features Ai) = 1) the probability is always one for any value of e. 

In Ref. [26j the authors have used an expression similar to Eq. (JTJ) , but with e = 1 and, as in Ref. [22] ]. they have 
modeled the strength of the VF as a probability of the interaction between agents and the VF. Increasing values of 
this probability implies a decreasing value of the probability for agents interacting between each other and then the 
corresponding diffusion of traits values between agents can be stopped which is not a realistic or desirable effect. They 
have found that only monocultural states arc obtained and only by introducing a noise rate it is possible to drive the 
system to a multicultural final state. 

In our case, we are not interested to include the effects of random perturbation effects, but we instead introduce 
another parameter related with the confidence of the information belonging to the VF. As "confidence" we understand 
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FIG. 3: Four possible cases of interaction for a system with _F = 3 features. Shared features are indicated inside a dashed 
rectangle. The probability of interaction p(i,j) (or p(i,M) in Case A) is indicating below. In each case the trait in a% x will be 
deleted by coping trait cr 3l (or by trait gm.\ in Case A). The probability of coping/deleting trait 1 is given by a) p' — C, b) 
p' = C, c) p' = 1 - C and c) p' = 1. 



here the credibility granted by agents to the information possessed by the VF. It is included as an extra probability 
C for agent i to copy an entry directly from the VF or an entry from another agent j that belongs to the VF either. 
It is also included as an extra probability 1 — C when agent i, when coping a trait, deletes and information the VF 
possesses in the same feature. 

To clarify this important concept we show in Fig. [3] four situations of interaction which resume all the possible 
cases. In Case A it is described an interaction between agent i and the VF which has been set to (0,0,0) without lost 
of generality. None of the nominal features are shared and then the probability of interaction p(i, M) depends only 
from the value of the effective features e according to the expression in the figure and is always larger than zero. In 
the practical case when agent i copies, for example the first entry, then the VF will be copied with probability p' = C 
which characterize the confidence of the information possesses by the VF. In the next cases, the interaction occurs 
between agents i and j which only share one trait of three possibles. The probability of interaction is then given by 
p(i,j) = 1/3 in all these cases. In Case B the nominal feature that agent i selects to copy from agent j coincides 
with the value the VF has in the same feature. Then, as in Case A, the corresponding trait is copied with probability 
p' = C. In Case C, when coping, agent i will delete its trait which is equal to that possessed by the VF. Then, it is 
deleted with probability p' = 1 — C. Finally, in Case D the traits copied and deleted are not related with the VF and 
then they are copied with probability p' — 1 . Note that according with these rules of interaction, when the confidence 
of traits possessed by the VF are C = 1, these traits are always copied with probability p' = 1 and never deleted. 
Otherwise, if the confidence of the traits possessed by the VF is C = 0, this traits are never copied (p' — 0) and are 
always deleted with probability p' = 1. Then, starting from the initial condition described above, the system evolves 
by iterating the following steps: 

(1) Select at random an agent i on the lattice, which is the active element. 

(2) Select at random, with equal probability, an agent of interaction. It could be one of the four nearest neighbors 
or the VF. 

(3) Calculate the overlap l(i,s) where s = j for the neighbor or s = M. for the VF. If s = A4, the agent i and the 
VF interact with the extended probability p(i, M). If s = j and < < F, agents i and j interact with 
the nominal probability p(i,j). 

(4) In case of interaction between agent i and agent s, choose a position trait h at random such that <ii h ^ o$ h (or 
Ci h ^ CM h ) an d then set Oi h = <jj h (or Oi h = <JM h ) according to: 

(4.1) if s = M. then set <ji h = <7m h with probability p' = C, 

(4.2) if s = j and <j, Jh = <Jm k then set Oi h = Oj h with probability p' = C, 

(4.3) if s = j and <Ji h = <Jm h then set Oi h — aj h with probability p' = 1 — C. 

(4.4) if s = j and both <jj h ^ <Jm h an d o% h ^ &M h i then set Oi h — <jj h with probability p' = 1. 
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Then, the full probability of agent i copies a trait from agent j is given by 

V(i,j) = ^p(i,j)p' (3) 
5 

while the full probability the agent i copies a trait from the VF is given by 

V(i,M) = -=p(i,M)p' (4) 
5 

Before studying the effects of the VF in our model let us review the original Axelrod's model. The computational 
dynamics of this model ends when the system reaches an absorbing state characterized by either j) — or l(i,j) = F 
for all pairs of closed neighbors A class of absorbing state, given by q F different configurations is called the 

"monocultural" state which corresponds to the case where j) = F for all pairs of closed neighbors. In this case, all 
agents in the network share the same trait at each feature (<Tj h = Oj H for all and the dynamics ends. Another 

class of absorbing state is called "multicultural" state and consist of at least two (or more) homogeneous domains 
which agents have cultural traits completely different. This way, two agents belonging to two different domains have 
zero overlap. The multicultural state is reached when each agent in the lattice has full of null overlap with all its 
neighbors. In these cases, a domain is given by a set of contiguous sites with identical state vector. 

It has been shown that the system reaches monocultural or multicultural states in dependence of lower (q < q c ) 
or higher values (q > q c ) of the cultural diversity q [l(j. To characterize the transition it has been considered two 
different order parameters: the average fraction of different cultural domains or the average number of agents in the 
biggest domain < S max > normalized to the number of lattice elements. 

In our case, it will be shown how our model produces a complex pattern of social behavior with affinities and 
repulsions to the influence of the VF in dependence of the effective features e and the confidence C of the information. 

III. NUMERICAL RESULTS. 

Numerical simulations have been carried out in lattices with L 2 = 30 x 30 agents and F = 4 features each. Different 
absorbing states have been found and we report the average realization number of agent in the largest domain over 
50 different initial conditions. 

Figure [5] shows the calculations of < S max > as a function of q at the absorbing state. Each panel shows the result 
for a certain value of the ratio e/F and different values of the confidence C in increasing order. The result using 
the Axelrod model without VF is included for comparison with full square dots. In this case, the system reaches a 
monocultural state at q < q c » 18 and a multicultural state at q > q c . In panel a) we set e = 1.0. It can be seen 
that when the value of the confidence C is small (C = 0.05), the results are very close to those of Axelrod model. 
The monocultural states remain unchanged at q < q c but the multicultural state is less "robust" and higher values 
of < S max > are obtained. For increasing values of C, higher values of < S max > for q > q c are obtained, as seen 
with C = 0.10 (hence the number of different domains in the multicultural state are smaller) and finally, at C — 
0.15, the multicultural states vanish and the system remains in monocultural states for all values of q. Then, it can 
be concluded that the increasing value of the confidence induces an homogenization of the cultural information that 
the system has, even at those values of q where the system reaches multicultural states when there is no VF. Then, 
multicultural states are unstable for increasing values of the confidence C at this value of the effective trait e. 

In Fig. H]b), c) and d) it is calculated the < S max > for smaller values of e. It can be seen that multicultural states 
at q > q c are again obtained for low values of C but now higher values are needed for the confidence to produce a 
cultural homogenization (< S max w 1 >). This can be seen when comparing the results for C = 0.15 in Fig. 2] a) with 
e = 1 and in Fig. 0]b) with e = 0.5. In Fig. 2]c) it can be seen that at a very small value of e (e = 0.1) the system 
has a multicultural state for q > q c with small increases of < S max > for increasing C, and even at e = 0.01 (Fig. [4] 
d)) the transition from a monoculture to a multiculture at q « 18 is independent of the confidence for values between 
0.0 and 0.2. 

Nevertheless, care has to be taken when analyzing the information of the induced monoculture at q > q c when the 
VF is present (as the case C = 0.15 in Fig. [4] a) and C = 0.20 in Fig. H]b)). It is interesting to know whether the 
greatest domain in the absorbing state characterized by < S max > has, or not, the information possesses by the VF 
in the corresponding nominal features. This information is indicated in Fig. 01 where all the full dots show absorbing 
states where the corresponding greatest domain does not possess the information of the VF in any of its features. 
That is, the overlap between VF and the cultural state of the largest domain is zero. Only at those absorbing states 
indicated by full stars (C = 1.00 in all panels and C = 0.50 on Fig. H]d)) the corresponding biggest domain fully 
shares the information at nominal features of the VF. Then, it is obtained the interesting result that at low confidence 
values C, the culture homogenization induced by the VF results in a negation or cancellation by the agents of the 
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FIG. 4: Calculation of the normalized average number of agent in the greatest domain at an absorbing state as a function of q 
averaged over 50 realizations. The values of the ratio e/F are a) 0.25, b) 0.125, c) 0.025 and d) 0.0025. At each panels, different 
value of the confidence C is taken into account in increasing order. The case of the Axelrod's model without VF is included in 
each panel with full square dots. Full rhombus indicate that the corresponding absorbing states do not share the information 
possessed by the external field, while full stars indicate full coincidence of the absorbing state with the VF. 



lattice of the information possesses by an external media. Here we call this phenomenon negative publicity effect, and 
it represents the process occurring in society when a group or different groups of people gather together, physically 
or intellectually, against an external action they consider misconceived. 

In order to study in more detail how the information of the VF spreads (or not) into the society, we have defined 
the parameter p which gives the percent of the total amount of the information that agents share in their nominal 
features with the VF as a function of time. It is given by 

i=i t=i 

and it is shown in Fig. [5l The calculation has been done for different values of the ratio e/F and different values of 
the confidence C taking in consideration only one initial state. In straight lines it is considered q = 34 while in dotted 
lines q = 8. When the dynamics starts from random initial conditions, the information possessed by the VF is already 
present in the lattice and shared by some agents. This anisotropy gives rise to the strong increase of its percent at 
the beginning but if the confidence value is low, the information is avoided by agents when traits are copied and it 
percent decreases with time after some maximum is achieved. This can be seen for C equal or below 0.39, 0.40 and 
0.41 in Fig. [5] a), b) and c) respectively. It is necessary a higher value of C for drive the system to a monocultural 
state with all agents aligned with the VF information. On the other hand, if the confidence value is high enough, the 
percentage of the VF information increases continuously until the absorbing state is reached. This can be seen for 
C = 0.40, 0.41 and 0.42 in Fig. a), b) and c) respectively. At the same time, there is a sharp discontinuity between 
the regions of confidence values where the system remains monoculture and multiculture, as can be seen between C = 
0.39 and 0.40 in Fig. [5] a), between C = 0.40 and 0.41 in Fig. \E\b) and between C = 0.41 and 0.42 in Fig. [He). For 
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Time steps / (L 2 . F) 

FIG. 5: Percent of the VF information in the lattice as a function of time for different values of the ratio e/F and the confidence 
C. In straight lines is used q = 34 while in dotted lines it is used q = 8. 

high values of C the system reaches the monocultural state in only few time steps, while for intermediate values of 
C the value of p changes slowly, at least until 2QQQL 2 F time steps where the simulation was artificially aborted. At 
those values of C the system seems to be in a quasi-stationary state and no absorbing state was found in these regions. 
For enough low values of C, as those shown in Fig. |4j the system reaches the multicultural state and there is almost 
no VF information on the lattice, as seen for C = 0.10 in Fig. [5] in the three panels. 

Finally, it is obtained that there are not fundamental differences between the behavior of the system with q = 8 
and q — 34 for higher values of e, as seen in Fig. [5] a) and b). Otherwise, in Fig. [5]c) it is reported an important 
difference in the dynamics for q = 8 and 34 in the region of confidence values where the system remains multicultural. 
This behavior is connected with the strong decrease of < S max > in Fig. 0]c) and d) for C = 1.00 and for C = 0.50 
in Fig. [4]d). In these cases, the high value of the confidence quickly drives the system toward the absorbing state 
because of the strong assimilation of the VF information by those agents that already share with it, at the beginning, 
the nominal traits. Then, the system quickly breaks down into domains with null overlap. 

IV. CONCLUSIONS. 

It is developed here a new model for the inclusion of an external vector field in the Axelrod model at zero temperature 
to describe the effects of the mass media on a social system. The clever design of publicity which allows the mass 
media to have influence over the whole society was included as a non-zero extended probability of traits being copied. 
This important effect is related with a parameter e which can be interpreted as an extra effective feature or features 
the mass media could have with all agents in society, beyond the nominal features. This effective feature(s) is(are) 
used by the VF to reinforce the frequency of certain information already present on the society or to introduce a new 
one. 
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It is also included in the model a confidence value of the information possessed by the VF. It is modeled as a 
probability of copying/deleting VF information which represents the criteria a person (or a group of persons) has 
(have) about what the mass media is proposing to the society. For very low values of this confidence, the dynamics 
recovers the Axelrod model with no external field, but an increasing value produces an homogenization on the society 
which would be multicultural without the external influence. This cultural homogenization is lined up against the 
acting influence of the VF with a zero overlap with the VF information. We have called here negative publicity to 
this effect. It simulates the behavior of people in society who gathers together against the external information they 
estimate wrong or incorrect. For large values of the confidence, the system reaches a quasi-stationary state with only 
slow changes of the amount of VF information in the society and no absorbing state was found for the amount of 
time steps tested. At higher enough values of the confidence, it is obtained a completed homogenized lattice lined 
up with the VF information, a situation that could be the purpose of mass media, politic party, etc. These results 
are qualitatively in agreement with the intuitive and realistic idea that with enough bombardment of the mass media 
information, that is accepted as valid and trusted, there will be a strong induced culture homogenization in the society. 
In other words, when people assume this information as personal the cultural differences tend to disappear. 

We also studied the dependence with the effective probability and it was found that at small enough values, the 
system is independent of the confidence value reproducing the results obtained in the Axelrod model. For increasing 
value of the effective probability, lower values of the confidence are needed to achieve monocultural systems lined up 
with the VF information. 

It is important to mention that the confidence parameter e can be experimental measured and, therefore, it could 
advance the results expected on a society when designing publicity. 
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